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We prove some new properties of fidelity (transition prob- 
ability) and concurrence, the latter defined by a straightfor- 
ward extension of Wootters' notation. Choose a conjugation 
and consider the dependence of fidelity or of concurrence on 
conjugated pairs of density operator. These functions turn 
out to be concave or convex roofs. Optimal decompositions 
are constructed. Some applications to two- and tripartite sys- 
tems illustrate the general theorems. 
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I. INTRODUCTION 

In Physics antilinearity is well known from symmetries 
with time reversal operations Q|, from second quantiza- 
tion, and from representation theory of groups and al- 
gebras. Quantum information theory offers several new 
interesting applications of antilinearity. In the present 
paper we are concerned with one of them. Antilinear 
operators are intrinsically non-local: One cannot tensor 
them consistently with the identity operator. They do 
not share the privilege of linear operators to allow 
execution in one part of a bipartite system while " doing 
nothing" in the other one. It seems, therefore, quite nat- 
ural to use antilinear operators to describe or to estimate 
effects of entanglement. Indeed, Hill and Wootters in ||] 
and Wootters in jj] used a particular conjugation, the 
Hill-Wootters conjugation, in order to get an explicit ex- 
pression for the entanglement of formation for two qubits. 
Their papers are the very starting point for the present 
contribution. I tried to distil a general method out of 
their proofs, and to construct explicitly the relevant opti- 
mal decompositions. The entanglement of formation con- 
cept is due to Bennett et all |5|. Also a peculiar basis, the 
magic basis, with which one can define the Hill- Wootters 
conjugation, is already in that important paper. 

In the 2-qubit case the entanglement of forma- 
tion is a function of just one other quantity, called 
(pre) concurrence, and the same optimal decomposi- 
tion of a state into pure ones can be used to calculate its 
entanglement of formation and its concurrence. In this 
form the statement becomes wrong for general states of 
a bipartite system different from the 2x2 case. But for 
density operators of rank two similar results seem not out 
of range. 

However, concurrence seems to be an interesting quan- 
tity in its own: It can be defined in higher dimensional 
Hilbert spaces and with respect to any conjugation 8 by 
an explicit expression (section II) which will be called 0- 



concurrence. It is a convex function on the state space 
(section HI), and it is a roof (see section V). Optimal 
decompositions can be obtained (section IV) in a con- 
structive manner, adding some news even for the 2-qubit 
case. Generally, the length of an optimal decompositions 
will be the smallest power of two which exceeds the di- 
mension of the Hilbert space. 

The idea, pointing to the definition of 0-concurrences, 
can be extended to another interesting quantity, to the 
fidelity, the square root of the transition probability 
0-fidelity as defined in section II, turns out to be a con- 
cave roof. Optimal decompositions can be gained simi- 
larly. 

The main proofs are in sections HI and IV. Section V 
is devoted to the roof concept, an interesting tool if 
combined with convexity or concavity. The last section 
contains some applications, mainly of 0-concurrences. 
There are conjugations in multipartite systems such that 
a non-zero 0-concurrence indicates inseparability. It is 
illustrated for bipartite (example 1) and for the 3-qubit 
systems (example 3). In a 2 x n bipartite system there is 
the possibility to bound entanglement of formation from 
below by the aid of 0-concurrences (example 2). After 
extending the method slightly (theorem 5) to a larger 
class of antilinear operators, example 4 treats 0-fidelity 
and -concurrence on some 2-dimensional subspaces of the 
2-qubit system. Though the result is essentially known 
for the concurrence, |3j, it explains a part of the method. 

Now I shortly call attention to some notations and 
rules, connected with antilinearity, to prepare what fol- 
lows below. An antilinear operator "i? acting on an Hilbert 
space H satisfies by definition 

t} (aiV'i + a2V'2) = a^V'i + a2^2 

If ■0 is an eigenvector of with eigenvalue A, eip is an 
eigenvector with eigenvalue e~^X for all unimodular num- 
bers e. The fact that the eigenvalues of an antilinear op- 
erator fill some circles in the complex plain will be used 
in the estimations of section HI. The product of two an- 
tilinear operators becomes linear, the product of an of 
antilinear operator and a linear one remains antilinear. 
The adjoint (or Hermitian adjoint), , of an antilinear 
operator is determined by the relation 

for all ^j, ip £ n. Notice = The standard rule 

{AB^ — for linear operators remains valid if one or 

both operators are replaced by antilinear ones. In par- 
ticular, with a complex number a and antilinear § one 
gets {a-dy = i?^a* = ad^ , i. e. taking the adjoint is a 
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linear procedure for antilinear operators. It follows: The 
set of operators which are antilinearly Hermitian (antilin- 
early self-adjoint), z9 = is a linear space of dimension 
d{d+ l)/2 if dimTi ~ d. Indeed, i) is antilinearly Hermi- 
tian iff {ijj, -dip) is symmetric. With respect to a basis the 
condition restricts the off-diagonal entries only. Complex 
diagonal entries are allowed. 

One calls antilinearly unitary or simply antiunitary 
iS 'd'' — i?^^. Basic knowledge about antiunitary op- 
erators is due to Wigner [Q. A conjugation, Q, is an 
antiunitary satisfying — 1. Writing O — 8^^ — 8^ 
shows the hermiticity (self-adjointness) of conjugations. 
Well studied examples are time reversal operators ||^ for 
Bose particles and for quantum systems with total inte- 
ger angular momentum . 

A conjugation 8 distinguishes in 7i a real subspace, Ti.Q, 
consisting of all 8-invariant vectors, i. e. of all eigenvec- 
tors of 8 with eigenvalue 1. No real subspace in TC is 
properly larger than Tig. Due to Hermiticity, &ip = tp 
and Qip = If result in 

so that the scalar product becomes real if restricted to 
Ti.0. In other words, Tie is not only a real subspace, it 
is a real Hilbert subspace. On the other hand, 8 can be 
gained as complex conjugation in every basis contained 
in He ■ This establishes a one-to-one correspondence be- 
tween maximal real Hilbert subspaces and conjugations. 
In a 1-qubit space, i. e. dimTi = 2, a conjugation in- 
duces a reflection of the Bloch sphere at a certain plane 
through its center. Selecting the 1-2-plane, the plane 
perpendicular to the 3-axis, as invariant plane, the effect 
of the conjugation to the Hermitian operator 

g ^ ^{xqI + xiai + X2(T2 + xsas), (1) 

that is g i-^ g = QgQ, reads 

Q ^ -^{xoi + xiai + X2cr2 - xsas). (2) 

Given a conjugation and a state vector, ip, we shall 
consider the absolute value of the transition amplitude 
between tp and Qip, or, what is the same, the square root 
of the transition probability between them. The quan- 
tity in question, |(^/;, 8'!/')|, is well defined for pure states. 
The problem addressed in the paper is to extend it to all 
states in a canonical way. In other words, we look for 
functions on the state space which are completely deter- 
mined by their pure state behaviour. This can be done 
by relying on the convex nature of the set of all density 
operators (states) which reflects the process of perform- 
ing Gibbsian mixtures, i.e. of convex sums. There is 
one and only one largest convex function coinciding at 
pure states with \{tl;,Qtp)\, and, following Wootters, I 
call it Q -concurrence. And there is exactly one smallest 
concave function within all functions which are concave 
extensions from the chosen values for pure states to all 
density operators. That function I call Q -fidelity. 



II. FIDELITY AND CONCURRENCE 

Let g and u be two density operators in an Hilbert 
space H. Their transition probability is denoted by 
P{g,uj), their fidelity^ the square root of the transition 
probability, is called F{g,Lu). It holds 

/P(^ = F(^),c^)=tr(V^^>V^)l/2 (3) 

Let Ti" be an ancillary Hilbert space. For any two vec- 
tors, (p,ip G Ti.<S) H"", which reduce to g and w, 

g = Tra\ip){(p\, w = Tra|'0)('0|, 

the transition amplitude is bounded from above by the 
fidehty, |(</j, V')] < F{g,uj). Indeed, F{g,uj) is the least 
number which fulfills this condition. Equivalently, as 
= P, a. suitably chosen von Neumann measurement 
in an ancillary system can cause a transition g ^-^ uj with 
probability P{g, oj). A larger transition probability, how- 
ever, is not possible The joined concavity of the fi- 
delity can be seen from 

F{g, Lo) = inf i ( iv{Xg) + iy[X-^Lo) ) (4) 

where X runs through all positive and invertible opera- 
tors X. A proof for finite-dimensional Hilbert spaces is 
as follows: Abbreviate by a and b the traces over Xg and 
X~^ui respectively. From Q one knows F"^ < ab. But 
2^/ab < a + b, and the right hand side of (^) cannot be 
smaller than the left one. If the density operators are 
invertible then there is a unique positive solution X of 

V V V -l/Z/ 1/2 l/2\ -1/2 

XgX — Lu, X — g [g ujg )g 

With this solution we get a = b and a = F, and (||) is 
saturated. Now we use continuity to extend the proof to 
all pairs of density operators. See als o_|lO[ . 
It is useful to extend the equations (^JT^lQ), and simi- 
lar ones to all positive operators with finite trace. The 
simple scaling properties of P, F, and related quantities 
make this is an easy task. Of course, the physical inter- 
pretation of P as a probability is bound to normalized 
density operators. 

Let 8 be a conjugation in an Hilbert space Ti. and 
abbreviate g :— OgQ. It is evident from (^) that 

Feig):=Fig,g) (5) 

is concave in g, [0. (||) will be called Q -fidelity of g. 

In order to introduce the (pre)concurrence ||] and Q| 
we need the ordered singular numbers, Ai > A2 > . . ., of 
yfg\JZj, that is 

{ Ai > A2 > . . . } = spectrum of {^/gL^Jy/gY^'^ (6) 

Having in mind Wootters' explicit expression for the en- 
tanglement of formation it is tempting to define for any 
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two density operators (whether normahzed or not) the 
function 



C{g, uj) := max{0, Ai — Xj} 



(7) 



and to call it concurrence of g and lo. 

A useful relation can be obtained if the rank of guj does 

not exceed two. Adding P = to the cross terms in 

the two non-vanishing eigenvalues cancel. But the sum 

of the squared eigenvalues (^) is equal to the trace of guj. 

Hence 

Cig, ujf + F{q, uf = 2Tr{guj) if rank [quj] < 2 (8) 

Finally, given a conjugation 0, we call Q-congruence of 
g the concurrence between g and its conjugate g, 



Cs{g) ■■= C{g,g), g^OgG 



(9) 



In contrast to the higher dimensional cases it is not 
hard to get explicit expressions if dimH = 2. With g 
given by (Q) and a conjugation acting as in (|^) one obtains 



(10) 



The next issue is to prove: Fq is a concave and Cq is 
a convex roof for every conjugation Q in every finite di- 
mensional Hilbert space. For the time being the finite di- 
mensionality of the Hilbert space is essential due to some 
unexamined mathematical problems in the case of infinite 
dimensions. Thus, in all what follows, dimTY = d < oo. 



III. PROPERTIES OF G-FIDELITY AND 
e-CONCURRENCE 

In this section we derive some implications from and 
start proving of 

Theorem 1 : Let 9 be a conjugation. Then 



Ce(f?) =min^|(0fe|e|0A 
Fe(f?) =max^|(0fc|e|0, 



(11) 



where the min and max has to run through all ensembles 
{01 1 02, ■ ■ •} such that 



Q - 



(12) 



is valid. 

The proof of the theorem will terminate in the next 
section. Up to that point we consider (|ll|) as a defini- 
tion of its left-hand-sides, and we shall draw conclusions 
without using (H) and (^ of the preceding section. 

Consider first the case g = Clearly, every de- 

composition ( p^ is gained by ipk = cik^p with numbers Ofc 
satisfying kfcP = 1- Hence 



(13) 



A simple consequence of ( |llD is homogeneity. For positive 
reals 

Cei^ig) ^ ^iCeig), Fei^ig) = iiFe{g), Vm>0 (14) 

Being in finite dimension the minimum (maximum) in 
( [III ) will be attained by certain decompositions (|lj). 
They are called optimal decompositions. 
Choosing optimal decompositions for Ce{g) and Ce{i^), 
their union is a decomposition for Ce(f? + w), though not 
necessarily an optimal one. Hence CQ{g) + C'q(lu) is an 
upper bound for CQ{g -\- to). Similar reasoning can be 
done for the 0-fidelity. Thus 



Ce{g + uj)<Ce{g) + Ce{iu) 
Fe(g + Lu) > Fe{g) + Fe{Lu) 



(15) 



showing subadditivity of Q -concurrence and superaddi- 
tivity of Q- fidelity. Because of its homogeneity ( p^ we 
conclude: 

Ce is convex, Fq is concave. 

Now we can go a step further, again without using 
arguments from the preceding section. Let f2 be the state 
space, i. e. the convex set of normalized density operators. 
If g is in this set, a decomposition ([l^ ) can be rewritten 
as a convex combination 



I0fc)(</'fc| 
{4>k\4>k) 



(16) 



Assuming that our decomposition ( p^ ) is optimal for, say, 
the 0-concurrence, we can write 

We conclude as following 0. Let C be another convex 
function on Q coinciding with C at the pure states. Then 
we have 

C'{g) < Y^PkC {tiu) = 5]pfcCe(^fc) 

But for a an optimal decomposition which of the G- 
concurrence the right hand sides coincides with Ce(p). 
A similar proof is for Fq . It results 

Theorem 2 : Ce is the largest convex function and 
Fq is the smallest concave function on the state space 
coinciding with \ {ip\Q\'il')\ at the pure states. 

To show that the right hand sides of (^l]) coincide with 
the definitions used in section 2, optimal decompositions 
will be gained in the next section. 



IV. OPTIMAL DECOMPOSITIONS 

In building optimal decompositions for our 0-fidelity 
and 0-concurrence the properties of antilinear operators 
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play a decisive role. Fix a density operator g and define 
an antilinear operator § by 

7?EE^, V^ey^ (17) 

Because 0t = 8, ?9 is antilinearly Hermitian. Hence 

Substituting (p = proves all the expectation values of 
■d^ real and not negative. Thus, d"^ is a linear positive 
operator and the same is with VW. Let us abbreviate 
Q = QgQ, so that can be written Remark, 
just to see what is going on, how the eigenvalues of the 
positive square root of d"^ have been used in section II to 
express Fq and Cq. 

Our next aim is to prove the existence of a conjugation, 
Oo, depending on g, with which we can polar decompose 



(18) 



^ = eoV7?2 = Vi92eo, r - ^g~g^. 



Let A^, A > be an eigenvalue of and Ti.'^ the 
Hilbert subspace of the corresponding eigenvectors. With 
-0 also dijj belongs to Ti^. Define on Ti^ the action 
Qq^) := X^^-Qi). On Ti.'^ the operator Oq is a conjuga- 
tion which commutes with i9. If one eigenvalue of i?^ is 
zero, Oo should induce on Ti^ an arbitrarily chosen con- 
jugation. Now Ti is decomposed as a direct orthogonal 
sum of Hilbert spaces of the form Ti,^ and Oq is given 
as an operator on every one of them. But this defines 
©0 uniquely as a conjugation on H, and ([l^ ) is proved. 
Choosing in every Ti^ a Oo-invariant basis, we get a com- 
mon eigenbasis, {-01, V'2, • ■ •}, such that 



By the help of (H) and (|l]) it is straightforward to prove 
the following, essentially known identities 



fe=l 



{vA^Wk) = ^Aj 



(24) 



The remarkable deviation from most uses of Hadamard 
matrices is in the appearance of the phase factors pro- 
duced by the antilinearity of 'd. They provide sufficient 
flexibility in adjusting the expectation values of d. By 
varying the Cj in the second equation arbitrarily, the ab- 
solute values of the numbers (</?fe|i?|<Pfe) fill completely the 
following interval of real numbers: 



d 
3 = 1 



\j I > max{0, Ai 



d 

i=2 



(25) 



Proof: a) The sum of the Xj is an upper bound (triangle 
inequality) and it is reached with ej"^ = 1 for all j. The 
simplest choice is = 1 for all j. b) If the Ai is not 
smaller than the sum of the remaining lambdas, a lower 
x\ where x is the maximum absolute 
\3 -I- . . .. Hence we get the asserted 
lower bound. The bound is attained for = 1 and 



bound is je-, ^Ai 
value of €2 ^2 + €3 ^^^3 



= i for j > 1. c) It remains to prove: If the as- 
sumption of b) is not valid, the lower bound should be 
reachable. In this case 



1?20jj. = \k1pk, 0O'0fc = Tpk 



(19) 



with ordered eigenvalues Ai, > A2, > .... 

The vectors constituting an optimal decomposition will 
be obtained by the help of real Hadamard matrices. They 
can be inductively gained by 



A. 



1 1 

1 -1 



A. 



(20) 



for m ~ 2,4,8,.... Let us denote by a^i the matrix 
elements of These entries are either 1 or —1. They 
fulfill 



(21) 



k=\ 



The number m is adjusted to the dimension c? of 7i by 

(22) 



m = 2"+\ 2"<dimH<2"+i 



With an arbitrary selection of d unimodular numbers 
(phase factors), ei, £2, . . ., we define with a basis (19) the 
vectors 



(fik ='^akieiipt, k = l,2,...,m 



(23) 



5^ Aj > Ai > A2 - ^ Xj 



i=2 



J=3 



The first inequality is the assumption, the second fol- 
lows because otherwise Ai < A2 in contradiction to the 
assumed ordering of the X^ ■ We like to conclude the ex- 
istence of a representation 



^i = iEs^'^^i 



as then the lower bound zero can be reached: We have to 
prove the same assertion as above, but now the length of 
the sum is d — 1. Hence the proof is done if ( p5| ) is true 
for sums of length less than d. Starting with d = 2, the 
proof terminates by induction to the length of the sum 
to be estimated. 

Given Ai,A2,... we choose unimodular numbers 
ei, £2, . . . saturating respectively the upper bound or the 
lower bound of (|2^). With this choice the vectors ( [2^ ) 
are denoted by (p^{to refer to the upper bound) and by 
ip'i^ (to indicate the use of the lower bound) respectively. 
From the construction follows that the insertion of 



fc = 1, . 



. , 771 



(26) 
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into ( p^ ) estimates ( pH ) as follows: 

d 

Ce < max{0, Ai - ^ Aj}, Fq>^ Xj 

These inequalities must be equalities. For the proof we 
use an arbitrary decomposition 1 = J2\Xk){xk\ of the 
unity, insert <j>k — y/QXk into and convert the sum 

to be estimated by the help of (|19|) into 

d 

k j=i 

At first we estimate concurrence by choosing the phases 
of tpj such that becomes real and positive. We 

get 

d 

a=|Ai-5|, 6s^A,^((0,,V^,)'l 

J=2 k 

for the sum in question. If |5| is larger than Ai we al- 
ready obtained a = with (jy^ = 4>~ ■ In the other case 
|6| cannot exceed A2 + A3 + . . ., i.e. a — |6| < Cq. 
Concerning the 0-fidelity the Schwarz inequality will be 
applied to the positive Hermitian form (0, Re- 
specting ( p^ and ( pj| ) one gets 

Therefore a cannot be larger than the trace of v^. The 
latter is equal to .Fe and we arrive at a < .Fe- ^ 
We have not only proved theorem 1 but also 

Corollary 3 : Let dimH = d and < d < 2"+^. 
For every g there exist optimal decompositions for the 
& -concurrence the length of which does not exceed 2""*"^. 
The same is true for the Q -fidelity. 

Remarks: a) Can the bounds for the optimal length be- 
come more stringent for certain dimensions of dimTi. 
The construction above seems to deny it. But a proof 
is missing, b) If d = 4 = 2 x 2, then n — 2 and there 
are optimal decompositions of maximal length four as 
shown by Wootters. See also jl^ for the optimal length 
problem. 

V. ROOFS 

We now call attention to some peculiarities of convex or 
concave function on the state space which admit optimal 
decompositions. These functions are quite different from 
unitarily invariant ones like, for instance, von Neumann 
entropy. The latter do not at all discriminate between 
pure states, they just estimate how strongly a state is 
mixed. Roofs, as defined below, and in particular convex 
or concave ones, draw all their information from their 
values at pure states. They try to interpolate between 



those values as linearly as possible. Let us see how it is 
achieved by two simple examples. 

In two d imensio ns Q-f idelity a nd 0-concurrence are 
given by \/l — a nd y/x^ + x^ on the unit ball xf + 

+ < 1 , see ( ]lO| ) . The first one remains constant on 
the planes X3 = constant, the second one does so along 
the lines Xi = ci, X2 = C2. The intersections of a plane or 
of a straight line with the unit ball are not only convex: 
The intersections can be gained as the convex hulls of the 
pure states they contain. 

In turning to the general case we denote by f2 the con- 
vex set of all normalized density operators on a finite 
dimensional Hilbert space and by the set of its ex- 

tremal points, i. e. the set of pure density operators. 
A convex subset of will be called a convex leaf of f2 
iff 

no = convex huU of (f^o n n^''"') (27) 

Let G — G{g) be a function on and ^Iq a convex leaf 
of ri. G is called convexly linear (or, equivalently, affine 
or fiat) on f2o if for all probability vectors pi,p2, . . . and 
for all choices of pure states 

7ri,7r2,... G llonl^P"™ (28) 

G satisfies the relation 

GiY,p,7r,)^Y.PM''^) (29) 

It is not necessary to check condition ( p^ ) for all possible 
convex linear combinations in case G is either convex or 
concave: 

Lemma R-1: Let G be convex or concave. If 

g^^qjTTj, qk>0 (30) 

is a decomposition of g into pure density operators 
TTi , 7r2 , . . . , and if 

G{g) = J2lM^,) (31) 
is valid then G is convexly linear on the convex hull of 

TTl, 7r2, . . .. 

Proof: Assume G is convex. Given g, there is convexly 
linear function I satisfying G > ^ on il, and G{g) = l{g). 
Together with (|3^) we get 

lig) = Gig) = J2 iM^j) > E *^('^^) 

Because the right hand term is l{g), the > symbol must 
be an equality sign. But G > I now enforces lijTj) — 
G{TTj) for the pure states involved in (|3l|). By the help 
of this equalities we estimate G{uj), oj — Pitti -t- . . ., by 

l{uj)<G{uj)<Y,PkG{7rk)^l{uj) 

and the inequality must be an equality. (The first in- 
equality sign is due to ^ < G, the second due to the 
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convexity of G.) This proves the lemma for convex G. 
Because — G is convex if G is concave, the lemma re- 
mains true for concave functions. Another proof is in 
@- 

By definition, G is a roof if Q can be covered by convex 
leaves such that G is convexly linear on every leaf of the 
covering. The covering is said to be a convex covering 
belonging to or compatible with G. 

There is a simple geometric picture beyond. Assume a 
real number g = gijr) is given for every pure state tt. The 
idea is to think of a wall, made of straight lines starting 
from TT and terminating at g(7r)7r. The demand is, to 
cover the state space by a roof, founded upon the wall, 
which is as flat as possible. To satisfy the demand one 
joins every two points on the wall by a straight line, every 
three points by a triangle, and so on. If the dimension 
of the polyhedra becomes large enough, (dimTi)^ in our 
case, the set of polyhedra covers fi, (an application of 
Caratheodory's theorem), and we stop. To get a roof we 
have to select a one-fold covering of from our huge set 
of polyhedra: There should be a function uj G{u>) such 
that X = G(w) whenever xlo is contained in one of the 
polyhedra of the selected covering. If it occurs, G{to)uj is 
a convex combination of the g{'Kj)-Kj which generate the 
polyhedron. Taking the trace yields a representation 



G{lo) 



m < (dim?^)^ 



From the bewildering manifold of roofs we select the high- 
est (or the lowest) one: Given u; we look for a polyhe- 
dron containing xlu with the largest (or with the smallest) 
possible real number x. Let us call this number G^{uj) 
respectively G^{uj). There is such a polyedron if g is 
continuous, because then the set of all polyhedra based 
on a bounded number of edges is compact. 
Some generalities can be abstracted from the construc- 
tion above, see 0. They are summarized in the 
following lemma. 

Lemma R-2 : Let g — g{'K) be a real and continuous 
function on the set of pure states. 

a) There is exactly one convex roof G^ and exactly one 
concave roof G'^ on f2 which coincides on f^P"™ with g. 

b) G"*" is the smallest concave function and G~ is the 
largest convex function which coincides at the pure states 
with g. 

c) It is 

G+{q)= max^pjg(7rj) 



g-{b) 



where the variations have to run through all convex de- 
compositions of Q with pure states. 

Starting the discussion above from g{ii) = 
where tt = ^-nd is a conjugation, we arrive at 



G"*" = Fq and G~ = Cq. Within the pure states belong- 
ing to one of the optimal decompositions of the preceding 
section the values .g(7r) remain constant. Hence G+ and 
G^ are constant on the convex leaf they generate: 

Corollary 4 : The Q -concurrence (respectively the 
Q -fidelity) allows for a convex foliation such that Gq (re- 
spectively Fq ) is constant over every of its leaves. 

As an immediate consequence, g — > f{CQ{g)) and 
g — !■ f{FQ{g)) are roofs over for every function f{x) 
defined on the unit interval. In general the roofs so ob- 
tained cease to be convex or concave. But there are some 
rules guaranteeing convexity (concavity) in some cases. 
To preserve convexity it suffices that / is convex and in- 
creasing. Concavity is guaranteed with / concave and 
decreasing jisll . 

Examples are: Cq with real 1 < s is a convex roof, Fq 
with < s < 1 is a concave roof. An important convex 
and increasing function, used by Hill and Wootters in 
and Q to get an expression for the entanglement of 
formation, reads 



fnwix) =: s( 



-) + < 



1- 



where s{y) abbreviates — j/lny. Thus 

Q ^ fHw{Ce{g)) 



), (32) 



(33) 



is a convex roof for every conjugation in every Hilbert 
space. However, only if the Hilbert space is 4- 
dimensional, and Q the Hill- Wootters conjugation, ( |3^ ) 
is equal to the entanglement of formation. In bipartite 
2 X 2n systems (|3|) can only be a lower bound to the 
entanglement of formation for appropriately chosen 0. 
(See the next section.) 

The following statements copy facts known in 2-qubit 
systems to a more general frame. 

The maximum of -Fe is one, and Fq = 1 is the equation 
of a convex leaf for Fe by lemma R-2. If FQ{g) — 1 
then SgS ^ g hy (|), and g has a basis of 0-invariant 
eigenvectors. 

The minimum of Cq is zero. The set of all states with 
vanishing 0-concurrence is a convex leaf with respect to 
Gq. If g is 9-invariant, QgQ ~ g, then Ge{g) = if and 
only if no eigenvalue of g exceeds 1/2. 

The entanglement of formation vanishes, as known 
from ||5[ , exactly for separable, i. e. classically correlated 
states [16|, Separability in a 2-qubit-system can 

equally well be characterized by the vanishing of Gq, O 
the Hill-Woottcrs conjugation. Again, just for 2-qubits, 
Fq = 1 is the equation for the convex hull of the maxi- 
mally entangled pure states. 



VI. EXAMPLES 

This section considers some possible applications of 
the general theorems. By looking at examples we ask 
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whether 0-concurrences can be used to decide separabil- 
ity problems in bi- and multipartite systems. In a 2-qubit 
system a density operator is separable if and only if its 
concurrence vanishes. Could one suppose similar state- 
ments in a higher dimensional or in a multi-qubit system? 
Certainly not with just one functional. But with suffi- 
ciently many it can work. Before treating the examples 
we have to return to a further issue in antilinearity. 

All conjugations of an Hilbert space are unitarily 
equivalent. From and the definitions of fidelity and 
concurrence one gets 

Fe'ig) - FeiUgU^), Ce'{g) = CeiUgU^) 

with Q' = WQU and every unitary operator U. How- 
ever, in a bi- or multipartite system, 



(34) 



one considers two conjugations equivalent iff there is a 
local unitary U such that 8' = U^QU. Some of these 
equivalence classes consist of tensor products of antiuni- 
tary operators. 



(35) 



To obtain a conjugation, the square of each factor must 
be a multiple of the appropriate identity, for example 
6^ = CqI". According to Wigner there are only two pos- 
sibilities, Co = ±1. Therefore, a factor in ( p5|) is either 
a conjugation or it is an antiunitary satisfying 6^ — —1. 
The number of the latter cases must be even to obtain a 
conjugation by (^. 

For the purpose of the present paper an antiunitary 9 
satisfying 0^ = — 1 is called a skew conjugation. While 
skew conjugations are mostly discussed in connection 
with time reversal of fermions, we need them as building 
blocks for conjugations in multipartite quantum systems. 
A skew conjugation fulfills = —6'' and 







(36) 



All expectation values of a skew conjugation vanish. 
There is a consequence for vectors ip G Ti. which are 
separable with respect to the first factor in (^^, say 

= (j)"' ® Lp. If the first antiunitary, 0^, is a skew con- 
jugation, the expectation value ('0, Qij:) must vanish. In 
other words: Let 6 be a conjugation ( p5[ ) and assume its 
first factor is a skew conjugation. If (V', 6'0) is not zero, 
\ili){'4>\ cannot be 7i"-separable. 

A skew conjugation, 9, allows for a representation 



(37) 



1 5: i < n, with a certain basis, V'l: "02, ■ • ■, called a 9- 
basis. By (g7|) the Hilbert space decomposes into a direct 
sum of 2-dimensional, ^-invariant Hilbert subspaces. Of 
course, any basis of H can serve as a 6'-basis for a certain 
skew conjugation 9. 



In 1-qubit spaces there is, up to a phase, just one 
skew conjugation 9 that may be defined by |0) —^ 
|1) — > — i|0). (The imaginary unit in the definition is 
by convention.) On the state space it induces the well 
known spin flip. With that definition 9 iS) 9 is the Hill- 
Wootters conjugation of a 2-qubit space. 

Example 1: Consider in (^) a direct product H — 
'H°'®Ti}' of two even-dimensional Hilbert spaces. We dis- 
tinguish a special class, JF, of conjugations: € if the 
conjugation can be written as the product Q — 9'^'' ® 9^ 
of two skew conjugations. Notice that, up to a phase, T 
consists of one conjugation in the 2-qubit case, the Hill- 
Wootters one. 

We have already seen from (36) that for this class 
(-0161-0) = if 0; is a product vector. Thus Cq{-k) = 
for every pure product state and for every Q ^ J- . But, 
as seen at the end of the preceding section, the equation 
C'e(p) = defines a convex leaf, i. e. Cq vanishes for all 
separable density operators. One may rephrase the state- 
ment by saying: If g is a state in a bipartite system and 
if we can find Q d J- such that Cq > 0, then g cannot be 
separable. 

We now complement the last statement: Let tt be pure. 
If C'q{tt) = is true for all Q Cz J- then tt is a product 
state, i.e. separable. 

For the proof we consider an arbitrary unit vector G 7i 
and assume dimW = 2n < dimTi''. We use the Schmidt 
decomposition 







ttl > Ct2 > 



(38) 



to define a skew conjugations in the two parts of our bi- 
partite system. 9a is defined by requiring 0°,027 ■ ■ ■ to 
be a 0a-basis. In TC'^ we complete, if necessary, the 0^ 
vectors to a basis which then is used as a the defining 9b- 
basis. After these preparations we consider Q = 9a ^ 9i,, 
a conjugation tailored for the vector (^8|). A straightfor- 
ward calculation yields 



leiV-) = 2^a2ja2,- 



(39) 



The sum on the right-hand-side can vanish only if all the 
Schmidt coefficients aj vanish with the exception of the 
largest one. Hence -0 must be a product state. 

Can we skip in the last statement the purity require- 
ment? It seems unlikely with the exception of the 2-qubit 
case. Thus we are faced with the problem to character- 
ize the set of states with vanishing 0-concurrences for 
all conjugations from JF. Let us call the set of all these 
states 57"^. As an intersection of convex leaves it is con- 
vex, but not necessarily a leaf. It contains all separable 
states. Moreover, a pure state is in if and only if it 
is separable. But not all extremal points of ^l'^ might be 
pure and, then, it will contain density operators which 
are not separable. 
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Example 2: We proceed with the setting of example 
1 and require W to be 2-dimensional. The requirement 
allows to bound the entanglement of formation from be- 
low for any even dimensional second factor TC' in the 
bipartite system. To do so we use (py) to establish 



2Vdetp = sup|(i/;,eV')|, 6 G J^, 
e 



(40) 



Here p denotes the partial trace of \ip){il))\ over the sec- 
ond factor, Ti}' . It then follows a lower bound for the 
entanglement of formation, E{q). 

E{g) > sup/w(Ce(e)) = /w(sup Ce((?)) (41) 
e e 

fnw is explained by (^). The equality sign is due to 
the monotonicity of Jhw- The right hand side of ( pl| ) 
is convex as a sup of convex functions of type (^3[). For 
pure states it coincides by (^0|) with the entanglement of 
formation. But the entanglement of formation is a con- 
vex roof by its definition, see |^ and point c) of lemma 
R-2. Hence the left hand side is the largest possible con- 
vex function with the described values for pure states. 

Example 3: Now we try a similar procedure as in 
example 1 for a 3-qubit-system. As already mentioned 
there is, after fixing a phase, only one skew conjugation, 
say 9, in a 2-dimensional Hilbert space. Every conjuga- 
tion in dimension two is of the form U9 with unitary U. 
H in (HJ) is now the direct product of three 2-dimensional 
Hilbert spaces. Consider the conjugations 



U9i 



U9®9, 9®9®U9 



(42) 



Let i/j £ H and G from this set. Then (^|8|^) is zero 
if 7/; is a product vector. A separable g allows for a con- 
vex decomposition with product states by definition. For 
Ce = determines a convex leave, Ce(^?) has to vanish. 
Turn now to the reverse and let be tt a pure states with 
Cq{tt) = for some conjugations listed in (|2|). The 
manifold of pure product states is 8-dimensional. We 
shall prove that eight equations Cq = with conjuga- 
tions from ( p^ ) are sufficient to decide whether is a 
product vector or not. 

This goes as follows. Write V as a sum |0)|</5o) + 

and start by the first set of conjugations listed in (42). 

We have to solve the equations 



The tilde abbreviates the Hill-Wootters conjugation 9(E)9. 
With unitaries of the form U\j) — ej\j) we see that ipk 
is orthogonal to ipk- Hence ipk is a product vector. To 
come to this conclusion, we need two diagonal unitaries. 
Next, with U equal to either cti or (T2, we see that ipo 
is orthogonal to ipi. Because both are product vectors, 
either the first or the second one of their constituents has 
to be orthogonal one to another. Hence, after checking 



Ce = with 4 conjugations from our list, we arrive, up 
to a local unitary, at one of two possibilities: 

|0)|</))|0) + |1)|</)')|1), |O)|O)|0) + |1)|1)|0') 

Choosing now a conjugation from the second group of 
(H) yields (<^|C/6'|0') = 0. We need just two of them to 
see that either </) = or 0' = has to take place, provided 
(j) is located at the second position in the direct product. 
To cover also the case with in the third position, we 
need two conjugations from the third group. 
Let TT be a pure state of a 3-quhit system. There are 8 
conjugations of the form such that n is a product 

state if and only if Ceiir) = is valid for all of them. 
It is tempting to ask whether one can prove similar state- 
ments for any multi-qubit system. I believe the answer 
is affirmative, but I did not check it. 

Last not least we are going to cure a curious shortcom- 
ing of the treatment in example 1: It cannot be applied 
if one of the factors of the bipartite system is odd dimen- 
sional: The set J- becomes empty. The same unsatisfac- 
tory event arises if no or only one factor of a multipartite 
system is even dimensional. 

Let us think, for example, the factor is 3-dimensional. 
To get an appropriate antilinear operator 9a we split W 
into a direct sum of a 2-dimensional and a 1-dimensional 
Hilbert space. In the former we equal 9a to a skew conju- 
gation. In the latter we set 9a to zero. We do not get an 
antiunitary operator, but an antilinear operator satisfy- 
ing 9^ = —9a. This relation suffices to guarantee (p6[). It 
seems natural, therefore, to allow in ( |35| ) the larger class 
of antilinear 9 fulfilling 9^ = ±9 as factors, and to require 
for the tensor product 8^ = 8 only. 
Returning to the bipartite system of example 1 we could 
consider the larger class of antilinear operators 
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9l 



so that 8 is antilinearly Hermitian, and, nevertheless, 
{tp, Qip) = for product vectors ^p. 

We arrive at the following general question: Do 8- 
fidelity §) and 8 -concurrence (O) remain concave respec- 
tively convex roofs for any antilinear self-adjoint 8. Go- 
ing through all the proofs one finds it essential that the 
antilinear operator := ^/qQ^/q is antilinearly Hermi- 
tian. For that reason one proves by literally the same 
arguments: 

Theorem 5 : Let 8 be antilinear and self-adjoint, 
8 = 8t. Then 

Fe :- F{g, 8^8), Cq := C(e, 8^8) 

is a concave respectively a convex roof. Theorem 1 and 
Corollaries 3 and 4 remain valid for them. 

Example 4: The final aim of the exercise is to deter- 
mine fidelity and concurrence of certain conjugated states 
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of rank two in in a 2-qubit space. The reader should con- 
sider the example as representative for a lot of others 
which need more calculation effort. 

Let 7^2 be a 2-dimensional Hilbert space. The transition 
probability can be given by elementary algebraic opera- 
tions [Q. For the present purpose an adequate expres- 
sion reads 



P{g,i^) = F{g,ijf =Trgiv + 2^/det~gdeUJ (43) 

By the aid of the equation can be converted to 

C{g, ujf = Trguj - 2^/AetgAeUl> (44) 

Let be an antilinear Hermitian operator acting on 7^2 . 
To get FQ{g) or Ce(g) we have to know the trace of 
gQgQ and the determinants of g and QgQ. 

After these preliminaries we think of TL2 as of a sub- 
space of a 2-qubit Hilbert space TL. We cannot use the 
Hill-Wootters conjugation Qhw in ( ^ ) or ( ^^ directly 
because, generally, 7^2 will not allow Qhw as a symme- 
try. Therefore we set 6 :— Q&hwQ with Q the projec- 
tion operator projecting H onto 7^2- ©, so defined, will 
be antilinearly Hermitian and it maps 7^2 into 7^2- By 
the little trick we see, abbreviating g = hw hw , 

FigrQ) ^ Feig), Cig,g)^Ceig) 

whenever g is supported by TL2 ■ 

We assume 7^2 is generated by two separable unit vec- 
tors. 



We choose their phases such that 



(45) 



(46) 



where a, b are positive real numbers between and 1. 
We get, by appropriately adjusting the free phase in the 
Hill-Wootters conjugation. 



(^i,ew^o) = {^o,eHwiJi) = v/(i-a')(i -fe') (47) 

We can replace Qhw by 9 = QQhwQ in (^ without 
changing its validity. Remind also that {ipi,Qipi) = 
because ■0^ is a product vector. 
We introduce a suitable basis by 



if' 



v/2(l + a6)^ 



^0 - ipi 
v/2(l - ab) 



By a short calculation one concludes from (47) 



a±ip 



a± 



V(l-a2)(l-62) 



l±ab 



(48) 



(49) 



Possessing a distinguished basis (48) in 7^2 we represent 



help of Pauli operators, see (|T]). The Pauli operators to 
the basis (^9|) are by convention 

a, = |^+)(^+| - \v^-){^-\,a, = + |(/.-)(^+| 

and (72 = io'io'3- Transforming g according to p — > QgQ 
can be accomplished by transforming the identity of 7^2 
and the just introduced Pauli operators. Using (p^. 



-0-3 



9 9 9 9 

al — at ai + a_ 



-1 + 



-as 



Q(7iQ = a_|_a_(Ti, 0(720 = — a+a_iT2 
One gets for the determinant 

det QgQ = (a+a_)^ det g 
and for the trace of gQgO 



al - a 



These expressions shall be inserted into (43) and 



Feigf = -[(a+ + a^)xQ + (a+ - a_) 



'2^3]^ — a_|_a_X2, 



Ce(£')^ = 4[("+ ^ "-)2^o + (a+ + a-)a;3]^ + a+a^xl 

One should have in mind xq = 1 for normalized density 
operators. Then the last equation represents just Woot- 
ters concurrence, C{g), for density operators supported 
by 7^2 • Returning to the amplitudes (Hq), a and 6, results 
in a more convenient form 

(e) = A^{xo~abxsY -{l-a^b^)xl (50) 
Ce{g) = A^{xs~abxof + {l~aW)xl (51) 



V(l-a2)(l-b2) 
1 - a262 



(52) 



One easily determines to convex leaves for these roofs: 
For Fq we fix a;o = 1, CC2, and x^ and let xi vary. We ob- 
tain a straight line in j;-space which intersects the Bloch 
sphere of 7^2 exactly twice, corresponding to the two 2:3- 
valucs with which xi , X2 , X3 becomes a unit vector. Along 
the line the 0-fidelity remains constant. 
The same procedure, however with fixing xo,xi,xs and 
varying 2:2, produces the convex foliation for the O- 
concurrence — which, in our example, is the Hill, Woot- 
ters one. 



any density operator g supported by 7^2 as usual by the 
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